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摘 要 : 针对 目前 反 向 计算 模型 还 无 法 实现 对 建筑 室内 边界 对 流 换 热 量 进行 反 向 计算 这 一 制约 性 差距 ， 采 用 温度 贡献 
率 方 法 ， 将 边界 对 流 换 热量 与 室内 测 点 温度 之 间 表 示 成 因果 关系 的 温度 贡献 因子 佐 阵 ， 基 于 计算 流体 力学 ， 将 最 小 二 
乘 与 Tikhonov 正则 化 方法 相 结合 ， 建 立 依据 室内 数 个 测 点 的 离散 温度 求解 边界 对 流 换 热量 的 反问 题 数 学 模型 。 应 用 三 
维 通 风 空 腔 和 某 建 筑 内 一 间 办 公 室 进行 实验 验证 ， 模 型 求解 值 与 实测 值 的 均 方 根 差 均 小 于 80%， 结 果 表 明 反 向 计算 模 
型 可 以 准确 对 室内 边界 对 流 换 热量 进行 求解 。 
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Inverse method to determine wall boundary convective heat fluxes in indoor environments 
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Abstract: Thus far, no research has inversely solved for boundary convective heat fluxes in indoor environments. This study 


proposed an inverse method based on Tikhonov regularisation and least-squares optimisation using computational fluid 
dynamics (CFD) to determine the wall boundary convective heat fluxes in indoor environments. The inverse model applied the 
CRI to establish a matrix to describe the cause-effect relation between the cabin wall boundary heat flux and the exhibited 
discrete temperature at certain points. To test and evaluate the proposed inverse model, the surface heat fluxes on walls in a 
three dimensional ventilation cavity and an office in a building were inversely solved. The RMS errors of the model were less 
than 80%. This study finds that the developed inverse method can accurately and efficiently determine the wall convective heat 


fluxes in indoor environments. 
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方法 反 向 求解 室内 颗粒 污染 源 的 释放 源 强 ; 文献 [2] 采 用 概率 密 
度 分 布 反 向 求解 通风 环境 内 的 污染 源 位 置 ; 文献 [3] 研 究 在 不 用 
在 设计 室内 环境 时 ， 必 须 考虑 热 边界 条 件 、 送 风 位 置 、 送 送 风 方 式 下 室内 污染 源 释放 位 置 对 污染 物 浓度 分 布 的 影响 ; 文 
风速 度 以 及 送 风 温度 等 环境 参数 。 为 获得 舒适 的 热 环境 ， 热 边 。 献 [和 在 不 同 送 风 方 式 的 基础 上 ， 又 增加 不 同 送 风速 度 这 一 条 
界 条 件 的 合理 设置 是 关键 因素 。 确 定 热 边界 条 件 的 常规 方法 须 。 件 ， 研 究 送 风 方 式 与 送 风 速度 不 同 组 合 情 况 下 对 室内 甲醛 浓度 
采用 一 种 相当 耗 时 的 迭代 式 的 猜测 与 校正 过 程 。 如 果 指 定 环 境 分布 的 影响 ; 文献 [5] 应 用 Fluent 数值 模拟 软件 研究 严寒 地 区 室 
内 的 温度 分 布 ， 那 么 就 需要 及 时 确定 与 之 对 应 的 热 边 界 条 件 ， 内 自然 通风 和 机 械 通 风 两 者 之 间 的 通风 效果 ， 并 分 别 研究 不 同 
这 样 将 极 大 地 加 快 热 环境 设计 过 程 。 这 是 一 个 基于 特定 目标 温 。 送 风速 度 与 不 同 送 风 温度 对 室内 不 同 高 度 位 置 处 甲醛 浓度 的 影 
度 反 向 确定 热 边界 条 件 的 计算 问题 ， 也 就 是 基于 有 限 的 实测 结 ” 响 。 
果 探 索 未 知 的 边界 条 件 ， 属 于 “由 果 及 因 ” 的 反问 题 范 畴 。 在 现代 工程 背景 下 的 传 热 反问 题 研究 有 了 一 定 进 展 ， 文 献 
许多 研究 人 员 在 室内 环境 中 进行 了 反 向 建 模 研究 ， 但 主要 ” [6] 构 建 一 种 反 演 方法 估算 烤箱 在 加 热 过 程 中 的 边界 温度 曲线 。 
集中 在 污染 源 辨识 0 问题 上 。 文 献 [] 采 用 正定 矩阵 因子 分 解 ”文献 [7] 提 出 一 种 基于 格林 函数 和 动态 观测 的 反 向 技术 , 在 钻井 
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过 程 中 对 工具 界面 的 温度 和 热流 量 进行 预测 ， 除 此 之 外 ， 由 于 
t 绒 梯度 法 的 计算 精度 较 高 外， 因此 共 轿 梯度 法 是 用 来 求解 传 
热 反 问题 的 一 种 常用 方法 。 文 献 [9,10] 采 用 共 力 梯度 法 反 向 求解 
混合 送 风 方 式 下 室内 边界 温度 变化 曲线 特征 ;文献 [11] 研 究 固 
体 热传导 反问 题 ， 应 用 共 轿 梯度 法 反 向 求解 固体 边界 导热 热流 
秆 ， 文 献 [12] 采 用 共 轿 梯度 法 反 演 高 炉 炉 衬 侵蚀 过 程 中 的 传 热 
线 特征 。 除 了 应 用 广泛 的 共 轿 梯度 法 ， 文 献 [13] 根 据 热 防护 
绝热 材料 氧化 铝 纤维 的 温度 测量 值 ， 采 用 遗传 算法 反 向 求解 氧 
化 铝 纤维 的 有 效 热 导 率 。 
上 述 反 问题 模型 均 可 利用 计算 流体 力学 (computational 
fluid dynamics, CFD ) 求解 ， 并 且 都 可 以 满足 计算 精度 需求 。 然 
而 , 这 些 方法 通常 代价 高 昂 , 需要 设 定 恰 当 的 迭代 循环 初始 值 ， 
通常 会 导致 计算 量 大 且 计 算 效率 低 。 而 且 上 述 的 反 模 型 均 基 于 
单一 的 热 边 界 条 件 进 行 反 向 求解 ， 没 有 考虑 当 多 个 热 边 界 条 件 
同时 存在 时 ， 各 个 热 边 界 条 件 之 间 是 如 何 影 响 温度 分 布 的 。 到 
目前 为 止 ， 还 没有 研究 针对 室内 环境 的 热 边 界 条 件 进 行 反 向 建 
模 ， 更 没有 一 种 便捷 高 效 的 方法 可 以 同时 求解 室内 多 个 边界 热 
流量 。 由 于 边界 对 流 换 热量 对 室内 热 湿 环 境 有 很 大 影响 , 因此 ， 
对 于 分 析 供 热 、 通 风 和 空调 系统 以 及 对 空调 环境 的 边界 传 热 特 
性 的 研究 具有 重要 意义 。 
反问 题 通常 都 属于 病态 问题 , 若 输入 数据 中 存在 微小 扰动 ， 
那么 会 对 整个 计算 过 程 产生 巨大 影响 。 常 采用 Tikhonov 正则 化 
方法 来 处 理 不 适 定 问题 以 提高 反 向 控制 方程 的 稳定 性 。 人 迄 今 ， 
于 室内 复杂 的 气流 模式 而 无 法 获得 能 量 方程 的 分 析 解 ， 因 此 
室内 热 环境 反 模 型 计算 结果 的 稳定 性 并 没有 完全 得 到 解决 。 除 
了 数学 优化 方法 ， 日 本 学 者 Kato04 在 1994 年 首次 提出 温度 贡 
献 率 (contribution ratio of indoor climate, CRI) 方法 ,该 方法 用 
来 分 析 室 内 单个 热 因 子 在 多 个 热 因 子 同 时 存在 时 如 何 影 响 其 它 
热 因 子 。 例 如 ， 室 内 墙 体 与 墙 体 之 间 的 辐射 换 热 ， 或 墙 体 与 室 
内 空气 之 间 的 对 流 换 热 。 然而 CRI 算法 只 考虑 对 流 换 热 对 室内 
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u-:gradT =div(T gradT)+ (1) 


其 中 : 7 为 空气 温度 ("C)， 也 为 空气 速度 (m/s), 厂 为 扩散 系 
数 (m/s)，Q 为 热源 (W)，Cp 为 定 压 比 热 (J/kg*C)，p 为 空气 
密度 (kg/m3)，div 为 向 量 散 度 ，grad 为 梯度 运算 。 

若 室 内 仅 存 在 一 个 热源 ， 那 么 由 该 热源 形成 的 稳 态 能 量 方 
程 如 下 : 


+ 

若 室 内 送 风 速度 恒定 不 变 ， 送 风 温度 维持 常 值 或 在 微小 范 
围 内 变化 ， 便 可 认为 室内 气流 模式 为 稳 态 流 场 ， 这 种 情况 下 ， 
仅 对 能 量 方程 进行 求解 即 可 得 到 温度 场 。 在 机 械 通风 条 件 下 ， 
强制 对 流 起 主导 作用 ， 若 热 边界 条 件 改变 时 引起 的 热 浮生 力 对 
流 场 造成 的 干扰 很 小 ， 可 假定 温度 场 是 一 个 线性 系统 nH， 室 内 
n(n>1) 个 热源 共同 作用 下 的 温度 场 等 于 单个 热源 形成 的 子 温 
度 场 的 线性 登 加 ， 表 达 式 如 下 : 


u:grad7 =div( 六 grad7 ) 


(2) 
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T(x90) T(x)+ By) ttD (2) -PT (ye) G) 


其 中 , 了 为 n 个 热源 下 的 总 温度 场 ，T; 为 第 i 个 热源 Ci 单独 存 
在 时 构成 的 子 温度 场 。 

由 于 室内 汝 流 特性 ,温度 场 在 严格 意义 上 不 完全 是 线性 的 ， 
但 与 平均 流 场 下 的 对 流 相 比 , 时 均 速 度 场 的 非 线性 显得 很 微小 ， 
因此 ， 假 设 室内 温度 场 的 线性 性 是 合理 的 。 
式 (3) 转 换 成 如 下 形式 : 


NY 


u':grad7 = div(T grad7)+ > 汪 (4) 


i=l Cp 


1.2 目标 函数 的 建立 
室内 热 边 界 条 件 的 改变 会 影响 温度 场 的 改变 ， 固 体 壁面 2 
间 仅 存在 辐射 换 热 ， 作 为 透明 介质 的 空气 不 吸收 固体 壁面 
的 辐射 射线 ， 因 此 ， 空 气 与 固体 壁面 之 间 的 对 流 换 热 是 室内 空 
度 


空气 温度 分 布 的 影响 ， 并 没有 考虑 墙 体 之 间 的 辐射 换 热 对 室内 

气 的 影响 。 因 为 辐射 换 热 只 在 国体 壁面 之 间 进 行 热量 传递 
于 空气 温度 ， 墙 体 壁 面 的 对 流 换 热 才 是 影响 室内 
空气 温度 的 最 重要 因素 ， 辐 射 换 热 仅 是 间接 因素 。 目 前 CRI 入 
法 仅 应 用 在 正 向 求解 室内 温度 分 布 的 问题 研究 中 ， 还 没有 将 访 
到 室内 热 环 境 的 反 向 求解 中 。 

综 上 所 述 ， 本 文 的 研究 目标 是 在 室内 环境 中 基于 CFD 和 
Tikhonov 正则 化 原理 ， 利 用 有 限 的 温度 测量 信息 反 向 求解 室内 
边界 对 流 换 热 量 ， 并 应 用 CRI 算法 加 速 反 向 计算 过 程 ， 最 终 根 
据 测 点 温度 快速 确定 室内 每 个 墙 体 壁 面 的 对 流 换 热量 。 


1 ”边界 对 流 换 热量 的 反 向 计算 模型 
1.1 温度 场 的 线性 系统 


气温 度 改变 的 主要 因素 。 根 据 流 场 稳 态 条 件 下 的 温度 场 满足 线 
性 性 ， 室 内 测 点 温度 变化 量 与 边界 对 流 换 热 变化 量 之 间 存 在 如 
下 关系 : 


T=AQ (5) 
其 中 : T=[T1，72，…，7T4] "为 测 点 温度 ，Q=[Q1,，Q2,…，Qn]™ 
为 边界 对 流 换 热量 。 由 于 假定 室内 流 场 回 定 ， 温 度 变 化 量 与 边 


界 对 流 换 热 变化 量 之 间 呈 线性 关系 , 因此 A 为 描述 边界 对 流 换 
热量 与 测 点 温度 之 间 因 果 关 系 的 线性 矩阵 。 

反 模 型 的 主要 目的 是 根据 室内 测 点 温度 求解 热 边界 条 件 ， 
于 矩阵 A 是 一 个 病态 矩阵 , 因此 反 向 求解 方程 (3) 是 一 个 不 
备 数值 稳定 性 的 病态 问题 。 鉴 于 此 ， 若 直接 对 方程 (53) 进 行 反 向 
求解 是 无 法 实现 的 ， 为 得 到 更 贴近 实际 情况 的 计算 值 ， 需 用 
定 的 方法 对 方程 (3) 进 行 特殊 处 理 ， 本 文采 用 Tikhonov 正则 化 


建筑 室内 的 热 边界 条 件 在 能 量 方程 中 为 热源 项 ， 
定律 推导 得 出 的 稳 态 能 量 方程 如 下 051: 


热力 当 


刘 
了 


方法 在 求解 过 程 中 增加 正则 化 项 来 加 强 计算 结果 的 稳定 性 。 通 
过 Tikhonov 正则 化 方法 将 方程 (5) 转 化 为 最 小 值 目标 函数 问题 ， 
表达 式 如 下 [171; 


Ew 
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Min 2(AQ,) =|A,,AQ, — AT + 8° LAQ,lb 


区 为 正则 化 矩阵 ， 


其 中 ，| 几 为 甜 阵 三 范 数 ，2 为 正则 化 参数 ， 


最 常用 的 表达 式 为 
1.2100. 
lo0o1210 a 
-oo 1 2 1 07) 
经 过 Tikhonov 正则 化 处 理 后 的 反 模 型 计算 4 者 果 如 下 : 
Q = (ATA + ALITL) x(ATAT) (8) 
根据 式 (8) 可 知 ， 若 已 知 测 点 温度 TT、 线 性 矩阵 A、 正 则 化 


矩阵 工 和 正则 化 参数 +4， 便 可 获得 边界 对 流 换 热 量 Q。 

如 何 选 取 合 适 的 正则 化 参数 4 是 Tikhonov 正则 化 方法 求解 
的 关键 ， 定 点 方法 (Fixed-point method，FP) 是 一 种 以 较 低 计 
算 代价 便 可 得 到 参数 4 的 简易 算法 [22-24], FP 算法 的 表达 式 如 
下 : 


,WR 0 9) 

根据 正则 化 修正 之 后 的 问题 (6)， 令 
x(4)=|AT -AQ (10) 
y2)=|z0ib (GD 


其 中 : x(0) 是 关于 4 的 单调 增 函 数 ,y(0) 是 关于 4 的 单调 减 函数 ， 
根据 这 一 性 质 , 可 以 得 到 如 下 关于 正则 化 参数 4 的 函数 表达 式 : 

开 (2)=x4)y (4) (12) 
其 中 : /是 参 变量 ， 由 于 函数 x) 和 y() 的 单调 性 ， 随 着 4 的 不 
断 途 代 ， 亚 (存在 极 值 点 ， 并 且 令 其 取得 极 值 的 4 就 是 优选 的 
局 部 正则 化 参数 。 对 亚 () 求 一 阶 导数 ， 令 其 等 于 零 ， 得 到 关于 
正则 化 参数 4 的 迭代 公式 入 ,=%( 和 )， 其 中 : 


办 (1)= sp 


[AT — 42 外 
lieil, 


对 于 变量 yj， 大 量 的 算 例 验 证 发 现 当 y=1 时 往往 可 以 得 到 
合理 的 正则 化 参数 4。 


(13) 


|A7-40D 
Zeal 


急 () 是 关于 4 的 单调 增 函 数 ， 它 具有 在 4 的 迭代 区 间 上 存 
在 一 个 是 点 和 一 个 四 点 的 几何 特征 ， 凸 点 所 对 应 的 4 就 是 局 部 
最 优 正则 化 参数 。 
1.3 温度 贡献 因子 A 的 构建 

采用 CRI 方法 获取 线性 矩阵 A 的 表达 式 。 热 边界 ; 对 测 点 
位 置 X 处 的 CRI 定义 式 03 如 下 : 


AT(X) _ 
AT 了 


用 (4)= 


(14) 


A7(X) 
[so/(pcw)] 


CRI,(X)= (15) 


其 中 : AT:(X) 是 第 i 个 边界 对 流 换 热 变 化 值 A 8; 在 测 点 义 处 引 
起 的 温度 变化 值 (*C)，ATP; 是 所 有 热 边 界 同时 存在 时 第 i 个 热 
边界 Qi 引起 的 温度 变化 值 (*C)，P 是 空气 密度 (kg/m3), Cy 是 定 
压 比 热 (J/kg*C), V 是 空气 通 量 (m3/s)。 
将 室内 单个 热 边界 独立 存在 时 通过 


对 流 换 热量 引起 测 点 温 


拆 


度 的 改变 量 记 为 ATs; 所 


闭 加 性 ， 


与 测 点 温度 变化 量 之 间 的 关系 式 如 下 : 


AT, 
AT, " Ca Cy 


的 表达 式 如 下 : 


Cy CD 和 
Co 0 


Ca Co 人 


将 表达 式 (15) 代 入 到 式 (16) 中 ， 得 到 温 


Ch Xi 


Cn) 0 


边界 对 
0 ... 0 
0 

0 pe X,, 


A=CX= 


2.1 实验 台 描 述 


9 热 边界 同时 存在 时 对 测 点 温度 的 贡 
献 率 记 为 矩阵 C， 单 个 热 边界 独立 存在 时 对 测 点 温度 的 贡献 率 
记 为 矩阵 入， 根据 温度 场 的 线性 


流 换 热 变化 量 
AQ, 
0 
AQ, 


度 贡献 因子 矩阵 A 


为 评价 反 


为 0.024 m 的 排 风 


ATCOY ADCAY ALC,pV 1 
AQ, AC, AO, CoPY 
ATC,OY ATCPVY ATCOY 6 1 1 
a | 
ATC,pV AT.C,pV ATC,pV ， ， 1 
AO A0, ~ Ao, CpV 
2 三 维 通 风 空 腔 的 实验 验证 
反 模 型 的 有 效 性 ， 搭 建 三 维 通风 空 腔 实验 台 ， 几 何 
尺寸 为 1.04 mx1.04mx0.7m(XxYx2Z)， 如 图 1 所 示 。 空 气 
左上 方 宽度 为 0.018 m 的 条 形 风 口 送 入 ， 从 空 腔 右 下 方 宽度 
排出 。 根据 图 1(b) 中 的 Pi-P4 共 4 个 测 点 温 
反 模 型 求解 图 1(a) 中 天 花 板 、 地 板 、 左 侧 壁 面 及 右 侧 


应 | 


谨 
Le 


壁面 的 对 流 换 热量 。 


N 
-rr---------- 
”~ 


大 化 板 


右 左 全 壁 抽 


(a) 三 维 空 腔 立 体 图 
1=1.04m 
DO0ISm Ee 
送 风 1 下 请 
一 
i 人 
Lo0.37 Mis 
VA) Mis 
T=15 °C 加 
A=0.00125m2s? P=200°C 可 
EE Ps=18.9°C T 
Py-187°C  w 
Ps=19.1°C 
Ei 


再 


| 图 


(b) 三 维 空 腔 中 心 截 


图 1 三 维 空 腔 通 


风 实 验 台 : (a) 立体 


图 ; (b) 中 心 截面 图 


)024 m 
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分 布 、 温 度 分布 以 及 污染 物 扩 散 的 有 力 工具 。 利 
数值 模拟 软件 对 三 
为 室内 气流 模式 为 消 流 ， 所 以 选取 RNG k-e 消 流 模型 ， 动 量 、 


.2 ”实验 验证 结果 


应 用 反 模 型 计算 的 前 提 条 件 是 稳 态 流 场 、 测 点 
I 量 值 、 热 物理 特性 和 空间 的 几何 尺寸 。CFD 是 模 


Pe 


通风 空 腔 的 速度 场 和 温度 场 进 


消 流 和 能 量 方程 的 离散 选取 二 阶 迎 风格 式 ， 采 月 


温度 、 边 界 
拟 室内 速度 


Fluent (6.3.26) 


行 模拟 。 
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3 ，” 某 建筑 办 公 室 的 实验 验证 


3.1 办 公 室 环境 描述 
为 验证 反 向 计算 模型 在 实际 应 用 中 的 可 行 性 ， 选 取 某 栋 建 
筑 内 的 一 间 办 公 室 进行 实验 测试 。 办 公 室 的 立体 图 如 图 2(a) 所 


日 Boussinesq 模 


示 ， 儿 何 尺寸 为 492 mx4.32mx2.42m (X xYx2)。 
办 公 室 的 一 个 角落 里 放置 一 个 文件 柜 ， 男 外 两 个 角落 里 放 


型 近似 热 浮 升力 。 


于 辐射 换 热 不 直接 作用 于 空气 


全 控制 在 1.0x10-3， 


在 模拟 过 程 中 不 考虑 辐射 换 热 模型 。 连 续 方程 和 动 


温度 ， 因 此 


置 两 张 桌 子 和 四 台电 脑 ， 桌 子 前 面 坐 着 四 个 人 ， 天 花 板 上 安装 


量 方程 的 收 


能 量 方程 的 收敛 控制 在 1.0x102。 在 CRI 
方法 计算 中 ， 存 在 4 个 热 边界 (天 花 板 、 地 板 、 左 侧 壁 面 、 右 侧 
壁面 )， 获 得 室内 稳 态 流 场 后 ， 便 可 获取 通风 空 腔 4 
的 CRI 分 布 。 根 据 4 个 壁面 的 CRI 分 布 ， 借 由 计算 式 (17) 构 建 
温度 贡献 因子 矩阵 A。 由 CRI 方法 构建 的 矩阵 描述 


个 固体 壁面 


了 壁面 对 流 


换 热量 与 测 点 温度 之 间 的 因果 关系 ， 由 于 该 矩阵 是 
类 型 的 函数 ， 因 此 不 受 稳 态 流 场 和 边界 热流 量 的 影 


P 


换 热量 。 应 用 

T 为 R4xl 维 向 量 ， 
正则 化 矩阵 。 前 已 述 及 ， 求 解 反 模 型 的 关键 是 要 
的 正则 化 参数 4。 一旦 确定 T、 工 和 14， 那么 根据 式 (8) 便 可 


AR 


中 


《全 
适 


4 个 测 点 温度 如 图 1(b) 所 示 , 在 z=0.35 m 截面 


1-P4 的 温度 值 为 已 


流 场 和 边界 
响 。 
中 心 线 上 的 


[信息 , 用 来 反 向 求解 4 个 固体 


壁面 的 对 流 


MATLAB 自 编 译 程序 对 反 模 型 进行 运算 。 


温度 贡献 因子 矩阵 AER™ 


求解 边界 对 流 换 热量 Q。 
Tikhonov 正则 化 方法 的 目的 是 将 修正 后 的 目标 函数 最 小 


化 。 应 月 


对 流 换 
知 ， 实 
0. 


对 流 换 热量 。 


， 工 采用 二 


找到 一 个 合 


将 4 带 入 到 式 (8) 中 便 


有 FP 方法 获得 4=1807.17 为 最 优 正则 化 参数 ， 如 图 2 所 
示 。FP 方法 计算 效率 高 ， 且 计算 能 力 强 。 
可 获得 4 个 壁面 的 对 流 换 热 量 ， 表 1 列 出 了 4 个 国 


体 壁面 实际 


据 量 与 反 模型 计算 结果 以 及 两 者 之 间 的 误差 
际 值 与 反 模 型 计算 值 较 接 近 ， 两 者 之 间 


。 从 表 1 可 


的 均 方 根 差 为 


6°C。 由 此 可 见 反 模型 可 以 确定 边界 对 流 换 热量 ， 且 具有 较 好 
的 实用 价值 。 最 终 根 据 测 点 温度 快速 确定 室内 每 个 墙 体 壁 面 的 


由 (4) 
3 
3 


2000 4000 6000 8000 10000 


1 


图 2 三 维 通风 


空 腔 反 模 型 


Ph Fixed-point (FP) 算 法 的 曲率 


表 1 三 维 通风 空 腔 边界 对 流 换 热量 的 实际 值 与 模型 计算 值 


热 边界 实际 值 (W) 模型 计算 值 (W) 误差 (W) 
天 花 板 -8.06 -7.53 0.53 
地 板 90.19 90.75 0.56 
左 侧 壁 -35.67 -36.34 -0.67 
右 侧 壁 -10.73 -11.39 -0.66 


主 :边界 对 流 换 热量 的 


FE 值 表示 热源 (释放 热量 ), 负 值 表示 热 


[( 吸 收 热量 ) 


六 瘟 日 光 灯 ， 地 板 上 安装 四 个 尺寸 为 0.5 mx0.5 m 的 散 流 器 ， 
天 花 板 上 安装 一 个 尺寸 为 0.3 mx0.3 m 的 排 风口 。 根 据 图 2(b) 
中 给 出 的 6 个 测 点 温度 ， 应 用 反 模 型 求解 图 2(a) 中 天 花 板 、 地 
板 、 东 墙 、 西 墙 、 南 墙 以 及 北 墙 的 对 流 换 热量 。 


国 


Pi=20.8 °C 
O 〇 


(b) 办 公 室 中 心 截面 图 
图 3 某 建筑 办 公 室 环 境 


1 (4) 
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图 4 ”办公室 反 模型 中 Fixed-point (FP) 算 法 的 曲率 
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3.2 ”实验 验证 结果 
采用 与 前 面 三维 通 风 空 腔 测试 实验 相同 的 CFD 模拟 方法 ， 
6 个 测 点 构成 的 了 为 Rex 维 向 量 ， 温度 贡献 因子 矩阵 A ERS*6， 
工 采 用 二 阶 正则 化 矩阵 。 应 用 FP 方法 获得 4=26.33 为 最 优 正则 
化 参数 ， 如 图 4 所 示 。 通 过 4 便 可 获得 6 个 墙 体 壁面 的 对 流 换 
热量 ， 表 2 列 出 了 6 个 墙 体 壁 面 实际 对 流 换 热 量 与 反 模型 计算 
结果 以 及 两 者 之 间 的 误差 ， 实 际 值 与 反 模 型 计算 值 之 间 的 均 方 
根 差 为 0.79*C。 说 明 反 模 型 在 实际 应 用 中 具有 可 行 性 。 
表 2 办 公 室 边 界 对 流 换 热量 的 实际 值 与 模型 计算 值 


热 边界 实际 值 (W) 模型 计算 值 (W) 误差 (W) 
天 花 板 -60.42 -61.23 0.81 
也 板 196.81 197.52 -0.71 
东 墙 76.41 77.19 -0.78 
西 墙 83.23 84.07 -0.84 
南 墙 55.09 55.85 -0.76 
北 墙 66.79 67.65 -0.86 


4 ”结束 语 


本 文 基于 室内 若干 个 测 点 温度 提出 一 种 确定 边界 对 流 换 热 


量 的 反 向 计算 模型 。 将 Tikhonov 正则 化 方法 与 最 小 二 乘 优化 结 
合 起 来 进行 反 向 建 模 。 正 则 化 参数 的 选取 对 反 模 型 求解 至 关 重 
要 。FP 方法 根据 过 代 函数 Bn(4) 的 局 部 凸 属性 构建 最 佳 正 则 化 
参数 。FP 方法 以 较 低 的 计算 成 本 得 到 可 靠 的 计算 结果 ， 这 是 
FP 算法 的 一 个 特点 , 这 使 得 实现 起 来 更 加 简单 和 高 效 。 应 用 三 
维 通风 空 腔 实 验 台 和 某 栋 建筑 内 的 一 间 办 公 室 对 反 模型 进行 验 
E， 实 测 值 与 模型 计算 值 之 间 的 均 方 根 差 均 小 于 80%， 结 果 表 
该 反 模型 能 够 准确 有 效 地 对 边界 对 流 换 热量 进行 反 向 确定 。 
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